Draft writing 



GENERATION OF A FIREBALL IN AGN HOT PLASMAS 

K. ASANO 

Division of Theoretical Astronomy, National Astronomical Observatory of Japan 

asanoOth . nao .ac.jp 
and 

F. Takahara 

Department of Earth and Space Science, Graduate School of Science, Osaka University 

takaharaSvega. ess . sci . osaka-u .ac.jp 
Submitted; accepted 

ABSTRACT 

Motivated by relativistic jets observed in active galactic nuclei (AGN), we sim- 
ulate outflows of electron-positron pairs strongly coupled with photons from nor- 
mal electron-proton plasmas. Using multi-fluid approximation and a Monte Carlo 
method of radiative transfer, we obtain spherically symmetric, steady solutions of 
radiation and pair outflows for the luminosity L < 10^^ erg s~^. For microphysics. 
Coulomb scattering, Compton scattering, bremsstrahlung, electron-positron pair 
annihilation and creation are taken into account. Although a significant amount 
of pairs outflow by powerful radiative force with a mildly relativistic velocity, 
the temperature is not high enough to avoid pair annihilation before the fireball 
becomes optically thin to scattering. Several caveats in the simulations are also 
discussed. 



Subject headings: plasmas — relativity — galaxies: jets 
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INTRODUCTION 



Relativistic jets are observed in active galactic nuclei (AGNs) and Galactic black hole 
candidates. The velocity of these jets is highly relativistic with a bulk Lorentz factor above 
10 and the kinetic power is almost comparable to the Eddington luminosity. The production 
and bulk acceleration mechanisms of these jets are still unknown, though many ideas have 
been proposed ranging from magneto-hydrodynamical to radiative ones. 

Although it is difficult to determine the matter c ontent of jets from observations, several 
independent argu ments favor electron-positron jets (ITakahara I 
19981 : iHoman fc Wardle Ill999l : iHirotani et al 



199d; 


Wardle et al. 


Il998l: 


Kino & Takahara 


booi 



199S 



2000; Hirotani 112005 



duced in accretion disks around the central black holes. Because the electron mass is much 
smaller than the proton mass, the produced electron-positron pairs can be more easily 
elected than protons. S ome papers discus s the accretion disk s with electron-positron outflows 



f lMisra fc Melia Ill995l : iLiang fc Li Ill995l : iLi fc Liang Ill996l : lYamasaki et al. 1119991 ). If the 



accretion disks form hot pair plasma strongly coupled with photons, the pl asma may be ther- 
mally accelerated like the fireball applied to gamma-ray bursts (GRBs) ( iRees fc Meszaros 
19921 ). As long as the initial conditions for t he fireball model are satisfied, the fiow acce l- 
erates to a relativistic velocity undoubtedly (jPiran et al. Ill993l : lAsano fc Iwamoto 1 12002| ) . 
However, the characteristic size of AGNs is too large to form fireballs which are in a complete 
thermal equilibrium at high temperatures comparable to the electron rest mass energy mgC^. 
In the standard fireball model, to make matters worse, electron-positron pairs are almost 
wholly annihilated in the course of the thermal expansion. 



To overcome this situation, llwamoto fc Takahara I (|2002| . 120041 ) showed that a "Wien 
fireball" , which is optically thick to Compton scattering but thin to absorption, results in a 
relativistic outfiow avoiding the difficulties of pair annihilation. When high-energy photons 
are provided sufficiently, copious electron-positron pairs are produced. From the size and 
luminosity of AGNs, the pair plasma is not expected to be in a complete thermal equilib- 
rium. The pair plasma may be optically thin to absorption but thick to scattering. Since 
photons and pairs are coupled by lepton scattering, the plasma can be thermally accelerated 
expending its internal energy. In this CclSG, clS long as the temperature remains relativistic 
^cC^)i sufficient pair creation by high-energy photons compensates for pair annihila- 
tion. Since the cross section of pair annihilation and that of Compton scattering are the 
same order, the number of pairs is almost conserved outside the photosphere, where photons 
and pairs are decoupled. Therefore, if the temperature at the photosphere is relativistic, a 
sufficient amount of electron-positron pairs survive as a relativistic outfiow. Although super- 
Eddington luminosity in spherical symmetry is needed in the Wien fireball model, the real 
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luminosity for collimated jets can be sub-Eddington. 



Pairs are formed via photon photon coUisions in the accretion disk composed of normal 
plasma (electron-proton). Assura i ng that pairs are confined with protons, several authors 



(ILightman 



198a 



1982 



Kusunose 



Kusunose 



Svensson 1982; Svensson 1984; Zdziarski 1984; Kusunose fc Takahara 



19871 ) investigated pair equilibrium plasmas. As is typically shown in 
(jl987l ). the plasma temperature becomes ~ O.lmgC^ for plasmas of size ~ 10^*^ 
cm and luminosity > 10"^^ ergs s~^, because of the high cooling rate of dense pairs. On the 
other hand, the Wien fireball model requires a super Eddington luminosity > 10^^ erg s~^ 
in the spherical symmetry and a high temperature > O.SmeC^. The series of studies of pair 
equilibrium plasmas imply that the equilibrium temperature is too low for the Wien fireball 
model. 

However, c onsidering the effects o f pair escape from the disk by their own pressure or 
radiative force (lYamasaki et al. Ill999l ). there is a possibility that the density and tempera- 
ture of pairs are drastically different from the values obtained assuming static plasmas. In 
order to prove pair ejection from the disk, one may need to treat multi-plasma dynamics 
with radiation field. However, there has been no study of formation and ejection of pairs 
from the accretion disk taking into account radiative transfer consistently. 

In this paper we simulate outflows of electron-positron pairs from hot plasmas with 
radiative transfer in order to discuss whether or not AGNs can produce fireballs with suf- 
ficiently high temperatures. For simplicity, the multi-fluid approximation is adopted, and 
the plasma is assumed to be spherically symmetric. The microscopic physical processes we 
take into account are Coulomb scattering, Compton scattering, bremsstrahlung, electron- 
positron pair annihilation and creation in electron, positron, and proton plasmas. In §2 we 
explain our method, and the numerical results are shown in §3. §4 is devoted to summary 
and discussion. 



2. METHOD 

In this paper we numerically obtain spherically symmetric, steady solutions of radiation 
and pair outflows from hot plasmas undergoing steady external heating. First of all we 
consider a static plasma consisting of protons (p) and electrons (e). The proton number 
density is assumed to be 

np(i?) = noexp(-(i?/i?o)'), (1) 

where R is the radius from the center. The total heating rate of the plasma is given as a 
parameter L. Assuming that the heating rate is proportional to Hp, the proton heating rate 
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per unit volume is obtained as 



Ho{R) = :^^^exp{-{R/Roy). (2) 

Inside this plasma electron-positron pairs (e^) are produced via photon photon collisions. In 
this calculation we assume the plasmas are divided into three fluids: proton {p), background 
electron (e), and pair (e^) fluids. The background electrons are conflned with protons to 
maintain charge neutrality. Thus, the number density of the e-fluid, Ug, is the same as that of 
the 1^ fluid, Up. On the other hand, positrons may move independently of the |>-fluid. In order 
to maintain charge neutrality positrons of the density n_|_ in the comoving frame accompany 
electrons of the same density ri- — n+. Namely, electrons are artificially divided into two 
components; one belongs to the e-fiuid, and another belongs to the pair-fiuid. Even if the 
pair-fluid flows with velocity (3 = v/c, p and e- fluids are assumed to be static. Although the 
validity of this multi-fluid approximation is not always assured, we adopt this approximation 
for simplicity. 

The conservation laws of energy, momentum, and number of electron-positron pairs in 
the rest frame of p-fluid K are given by 

l,^[R:'(,.^p.}ru]^H, (3) 
i_^[«Vf/]=0. (5) 

where n± = 2n+, e± and P± are the number density, energy density and pressure of the pair- 
fluid in the fluid rest frame K', respectively. The velocity of the pair-fluid is represented 



by the Lorentz factor, F = 1/ yl — and 4- velocity U = \/r^ — 1. The pressure is 
written as n±T±, where T± is the temperature of the pair-fluid. Normalizing T± by the 
electron mass as 9± — T±/{meC^), the energy density in the rest frame, e±, is written as 
n±m^c^{K^{l/9±)/K2{l/9±) —9±), where is the ith order modifled Bessel function of the 
second kind. The source terms H, F, and Q, which are measured in the coordinate frame, 
are explained in the following subsections. 



2.1. Energy Budget 

The energy source term for pair- fluid H is expressed as H — Hp + + H^, where Hp 
and Hg are the heating rates by protons and background electrons via Coulomb scattering, 
respectively, H^ is the heating rate due to interactions with photons. 
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When there exists relative velocity betwe en the pair-fluid an d background fluids, it 
may be hard to obtain Hp and analytically. lAsano et al. I (120061 ) numerically calculated 



Hp and H^. for two Maxwell-Boltzmann gasses with relative velocities U = 10~^-10^ using 
relativistic formulae. Since these rates are functions of U, T±, and background proton 
(electron) temperature Tp (Tg), we have no simple fltting formulae for Hp and Hg. However, 
if [/ < 1, numerically obtained Hp agrees rather well with the analytical formula of t/ = 0, 



Hn 



+ 0. 



n±np- 



3me<JT 



2mp K2{l/e±)K2{l/ep[ 



In A X 



e± + 9p 

0+6„ 



2Ka 



d± + dp 
9±9p 



(6) 



( IStepney fc Guilbertlll983l ). where 9p = Tp/{mpC^) and In A is the Coulomb logarithm (we 
adopt In A = 20 here after) . In our s imula tion we use the above formula for Hp and tables of 
the case of f/ = in lAsano et al. I (120061 ) for Hg, irrespective of U. As will be shown later 
in our results, U becomes ~ 1 in low-density regions, where the Coulomb energy transfer is 
not effective. Therefore, we neglect the effect of flnite U on Hp and H^ in our case. 

The heating rate of pair-fluid through photons H-^ is divided into two parts as H^ = 
Heir — Hrad, where H^ad is the cooling rate due to photon production (bremsstrahlung and 
pair annihilation), and Hen is the heating rate by circumambient photons via Compton 
scattering and pair-creation. In this paper we only consider thermal bremsstrahlung and 
pa ir-annihilation fo r photon-production processes. We use the formulae of the cooling rates 



m iSvenssonI (119821 ) for bremsstrahlung {ep, ee, and e+e") and pair annihilation. The ra- 
diative cooling rates due to scattering of two particles in the background fluids {ep and ee- 
bremsstrahlung) are easily obtained. Hereafter, we call such emissions "background-fluids 
emissions". Given the temperature Tg, these cooling rates via bremsstrahlung, and B]^f 
[erg cm~'^ s~^], are proportional to n^, where the superscript "bg" designates the background 
fluids. 



On the other hand, bremsstrahlung emission and pair annihilation of two particles in 
the pair-fluid (pair- fluid emissions) are calculated in the comoving frame K'. The Lorentz 
transformation implies that the cooling rates due to bremsstrahlung, B^" , and B^"=B^" , 
and the rate due to pair annihilation, A^'^''' in the coordinate frame K are F times those in 
the frame K'. 

We should notice that there are emission processes between two particles of different 
fluids; one belongs to the pair-fluid and another belongs to the background fluids. We 
call such emissions "inter-fluid emissions". Since the pair-fluid is moving velocity, inter- 
fluid emissions are not isotropic in any frame. In our simulation we adopt the following 
approximation for the cooling rates due to inter-fluid emissions. In the frame K^,, which 
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moves with velocity = f//(l + F), the two fluids flow symmetrically; the background 
fluids flow with velocity while the pair-fluid flows towards the opposite direction with the 
same speed (3^,. In this frame the densities of the e-fluid and pair- fluid are np^, — Y^^Hp and 
77,-1-^, = r^,?T,-|-, respectively, where T^, = 1/ ^/l — We use the usual formulae for thermal 
and isotropic plasma adopting the temperature as max (Tav, /5*r^,) and densities n^* and n±*, 
where Tav = ^*{T± +Te)/2 for bremsstrahlung, i?g°* and B]^, and pair annihilation, A™*, 
or Tav = ^*T± for bremsstrahlungs, B^^^ = B^^^. The superscript "int" means inter- fluid 
emissions. The cooling rates in the frame K become F,, times those in the frame K^:. We 
assume that two fluids lose energies equally via these emissions, except for B^^^ and B^^ 
processes, in which only the pair-fluid loses energy. 

The above approximation may seem rough. However, as will be shown in our results, 
n± ^ Hp and T± ~ Tg are plausible for energetic jets. Therefore, inter-fluid emissions have 
minor effects on the motion of the pair-fluid. 

From the above assumptions we obtain 

Hrad = - f AP'" + + 25^^ + 25^7+ 

c \ ^ 



The heating rate Hdi by circumambient photons via pair creation and Compton scat- 
tering is obtained from calculation of radiative transfer, which is explained in the next 
subsection. 

The heating rate of background electrons. Hep, by protons is also obtained from ([6]) 
replacing T± and n± by Tg and Up, respectively. Then, we can obtain the temperature of 
protons Tp by energy balance equation as 

Hq = Hp + Hep. (8) 

The background electrons are also heated by photons via Compton scattering. We write the 
heating rate of e-fluid by circumambient photons as He^dv Then, the temperature of e-fluid, 
Te, is obtained from 

Hep -He+ He,cir = ^ Q^^"* + + 5^^ + 

^BT- + ^BT:) ■ (9) 
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2.2. Radiative Transfer 



We solve radiative transfer witli tlie Monte Carlo metliod, taking into a ccount pair- 

creati on and Compton scattering. Our numerical metliod is basically the same as llwamoto fc Takahara 
( 120041 ) ■ We shortly explain the method in this subsection. 



We divide the calculation regime into many spherical shell elements. From each shell, 
photons are emitted via background, pair, and inter-fiuid emissions, which are assumed to 
be isotropic in the frames, K, K', and K^,, respectively. For bremsstrahlung, we adopt an 



approximation of number spectra oc exp [— (2e/T)°-^], where T = Tg, T-t, or Tav for each 
emission process, normalized by the emissivities discussed in the former subsection. This 
approximation well reproduces the specral shapes o f brenasstrahlung numer ically obtained, 
based on the methods of quantum electrodynamics (lHaugill975l . Il985l . |2003| ). On the other 
hand, we adopt a fitting formula of IZdziarski I (jl980l ) for pair-annihilation spectrum. 



The pair-fluid emission and the inter-fluid emission are assumed to be isotropic in the 
frames K' and K^,, respectively. We define the emission coefficient, [cm~'^ s~^ sr^^]: 

emitted energy per unit time per unit volume per unit solid angle and per unit photon-energy, 
where is the cosine between the direction of the photon and the radial direction. Taking the 
inter-fiuid emissions as an example, th e emission coefficient in the frame K is obtained from 
the Lorentz invariant j^/e^ = je^/ei (IRybicki fc Lightman I Il979l ) and = eV^{l — P^ji), 
where is the photon-energy in the frame K^,. For the pair-fluid emissions too, we obtain 
the emission coefficients in the same way. 

Generating photons according to the above method, we simulate the trajectory of each 
photon propagating through the three fluids and photon field, where the photon interacts 
with electrons, positrons, and photons themselves. Elementary processes to be considered 
are Compto n scattering and pair creat ion. We use the Klein-Nishina formula and the repre- 
sentation in 



Gould fc Schreder I (119671 ) for the cross sections of Compton scattering and pair 



creation, respectively. We simulate all these processes by th e Monte Carlo method . Deta ils 
of the method of following the photon trajectory is given in llwamoto &: Takahara I (120041 ). 



After simulating all photons emitted from all shells toward all directions, we numerically 
obtain the photon distribution in the phase space, the heating rates i^cir and Hg cn, the 
radiative force on the pair-fluid F^, and the pair creation rate Q^. However, results depend 
on the photon distribution assumed in advance. Therefore, given the parameters of the three 
fluids, we simulate iteratively until the photon distribution converges. 
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2.3. Momentum and Number Sources 

The momentum source F is expressed a.s F = F^ — F^-ca — Fg — Fp, where F^ is the 
radiative force due to pair creation and Compton scattering, Frea is the reaction force due 
to photon-emission, and F^ and Fp are the frictional forces due to Coulomb scattering with 
background electrons and protons, respectively. As was mentioned in the former subsection, 
F^ is obtained from the calculation of radiative transfer. 



Asano et al. I (120061 ) calculated Fg and Fp numerically for two Maxwell-Boltzmann 
gasses with relative velocities U = 10~^-10^ using relativistic formulae, and they express 
their results as a form 

Fe^p = (Tn±)npmeC^fe,p(T±,Te,p,U), (10) 
where rn+ is the de nsity of pairs in the frame K. The functions f^^p are listed as tables in 



Asano et al. I (120061 ). For U <^ 1, f^^p increases with U monotonically. Beyond f/ ~ 1, 
decreases with f/, while /e is nearly constant. The lower temperatures T± and Te become, 
the larger the frictional forces become. 

Although photons are emitted isotropically in the frame K' , the photons are beamed in 
the frame K and the pair-fluid loses momentum. The Lorentz transformation implies 

Frea = ^ [A^^' + 2Blt + , (11) 

where we neglect the contribution of the inter-fiuid emission. In the strict sense the inter- 
fluid emission causes momentum transfer between the pair-fluid and the background fluids. 
In order to obtain the momentum transfer in bremsstrahlung, the scattering processes of two 
particles with photon emission, we need messy and heavy calculations using the differential 
cross sections as functions of both the particle scattering angle and the photon angle. On 
th e other hand, t he mo mentum transfers due to pair-annihilation are numerically calculated 
in lAsano et al. I (120061 ) . which shows that the rate due to pair-annihilation is negligible in 
comparison with that due to Coulomb scattering. The momentum transfers due to processes 
accompanying photon emission, including bremsstrahlung, may be negligible. The smallness 
of momentum transfer and the difficulties in calculations are reasons why we have neglected 
the contribution of the inter-fluid emission in equation (ITTl) . 

Of course, the background fluids are also pushed outward by photons and the frictional 
force. In this paper, we assume that the background fluids are maintained to be static; such 
forces are balanced with the gravitational force by the central black hole, which is implicit 
in our geometry of computation. The e-fluid suffers from the counteractive frictional force 
Ff. and the radiative force, which can be stronger than the total force on the ]9-fluid. There 
should be a certain degree of charge separation between the e-fluid and p-fluid to maintain 
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the e-fluid static. The resultant static electric field may affect the motion of the pair-fluid. 
Electrons in the pair-fluid are decelerated, while positrons are accelerated. The force due to 
the static electric fleld on the pair-fluid becomes zero on the average. As long as we adopt 
the multi-fluid approximation, we can neglect the effect of the static electric fleld, although 
there is possibility that the kinematical effects of the electric fleld heat the pair-fluid. 

The source term Q, which is Lorentz invariant, is expressed as Q = — Qanm where 
and Qann are the pair creation and annihilation rates, respecti vely. Whi l e i s obtained 



from the calculation of radiative transfer, we adopt a formula in ISvenssonI (119821 ) for Qann- 
As is the case of radiative cooling, the positron annihilation rate is divided into two parts as 
Qann = Qann + Qann- The method of calculatiou of Qann the Same way as the inter- fluid 
cooling, in which n^*, n±^,, and T^^ are used. Even if a positron annihilates with an electron 
in the e-fluid, the number density of e-fluid should be invariant to maintain the multi-fluid 
approximation. Therefore, we should consider that an electron in the pair-fluid is traded to 
the e-fluid in this case. As a result two particles in the pair fluid disappear in each inter-fluid 
annihilation. The trade of electrons means that energy and momentum are also transferred 
between the two fluids. However, as we have mentioned several times, such energy and 
momentum transfer may be negligible in this case too. Although this unavoidable treatment 
in the inter-fluid annihilation indicates a limit in the multi-fluid approximation, we adopt 
this approximation to save cost in computation. 



2.4. Flow Solution 

We solve n±, F, T-t, Tg, and Tp from equations ©-([SD,®, and First of all we set 
plausible values for the three fluids. Then, combined with radiative transfer to obtain the 
source terms, we mimic the time evolution of the plasmas. Since the photon diffusion time 
scale is much longer than the dynamical time scale, we cannot follow the time evolution 
of both the plasmas and photon fleld consistently. However, the flnal steady solutions are 
consistently obtained by mimicking the time evolution for a sufficiently long time. 



3. NUMERICAL RESULTS 

In our numerical simulation, the energy space of photons x = e/me(? is divided into 
70 bins in the logarithmic manner between 10~^ and 10^. For the direction of photons, /i is 
divided into 20 bins between —1 and 1. The outer boundary is set aX R = i?out = 2-Ro) and 
the radial coordinate is divided into 40 shells. We follow more than 10^ paths of photons to 
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Fig. 1. — Positron density (solid, left axis) in the frame K' and 4-velocity U (dashed, 
right axis). The proton density is also plotted for reference. 



solve radiative transfer. Some figures in this section are plotted using results after following 
5.6 X 10^ paths of photons. We do not take into account photons from outside of the boundary. 
However, such photons may not affect the outflow as will be shown later. 



In this calculation the parameters are the radius Rq, density uq and luminosity L. For 
several values of L above 10^"^ erg s~^, we simulate plasmas for two sets of parameters: 
Ro = 10^^ cm with uq = 10^° cm~^, and -Rq = 3 x 10^^ cm with no = | x 10^° cm"^. These 
values are characteristic of accretion plasmas in AGNs. In both cases, the pair free Thomson 
scattering depth Xp ~ no-RoC"T ~ 0.7. 



When L is less than 



]^q45 gj,g g 1^ ^YiQ results are almost the same as those in lKusunose 



(119871 ). As the luminosity L increases, the temperature of pairs decreases, while the pair 
density increases. The positron density is comparable to the proto n density. Howeve r, a 
depression of the pair density in the central region, which is seen in iKusunose I ( 119871 ) for 
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Fig. 2. — Temperatures 9± (solid, left axis) and Op = Tp/{meC^) (dashed, right axis). 



L ~ 10^^ erg s~^, is not observed in our simulation. The pair density decreases with radius 
R monotonically. The pair fluid escapes from the boundary by their own pressure. The 

radiative force is not effective to drive the outflow in these cases. However, the energy 
fraction of the outflow is negligible compared with the total luminosity L. The almost all 
energy injected in the plasmas escape as photon energy. Outside the boundary, most of these 
photons may escape without interaction with the pair-fluid. Therefore, the pair-fluid does 
not behave as a flreball in these cases, which is out of our interest. 

As the luminosity L increases above IC^^ erg s^^, the pair density overwhelms the proton 
density. Although the pair temperature 9^ continues to decrease with growth of L, Oj,- is 
always larger than 0.1. For L > 10'^^ erg s~^, a considerable amount of energy escapes in the 
form of the pair-fluid. In this high-luminosity case, the outflow is driven by radiative force 
rather than their own pressure. 
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In this paper we show detailed results for L = 10^^ erg s~^, i?o = 3 x 10^^ cm, and 
^0 = ^ >^ 10^° cm~^ as a representative case. The positron density and 4-velocity U are 
plotted in Figure 1. The ratio of number densities of positrons and protons z = Tn^/n^ 
is 74 and 440 at the center and the boundary R = Rout, respectively. The 4-velocity U 
monotonically increases with R, and becomes 1.01 at the boundary. The mean free path for 
the Thomson scattering at the boundary is It = T/{n±aT) = 1.1 x 10^^ cm, which is much 
shorter than R = i?out = 6 x 10^^ cm. Therefore, even outside the boundary, the outflow 
and photons are tightly coupled and behave as a fireball. In other words, our simulation 
is confined well inside the photosphere, which raises some problem with the treatment of 
the radiative transfer. However, as is seen later, since the photons are beamed in the outer 
region, the problem may not be so severe. 

The outflow may continue to be accelerated outside -Rout- As conventionally used in the 
GRB standard model, we define the ratio of the total luminosity to the mass ejection rate 
as 

^ 4:7iRl^^n±{Ront)U{Rovit)meC^' ^^'^^ 

In the GRB models this value, which corresponds to the final Lorentz factor of the outflow, 
is defined using the baryon ejection rate. In our case too, if the number of pairs is practically 
conserved within the photosph eric radius, the final Lorentz factor of the flow is on the same 



order of f], as demonstrated in llwamoto &: Takahara I (120041 ). Our simulation shows ~ 34 



which is enough to explain the high Lorentz factor of AGN jets. However, in order to conserve 
the number of pairs, a relativistic photon temperature at the photosphere is required. 

The temperatures 9± and 6p = Tp/{mi,c^) are plotted in Figure 2. At the center the 
pair temperature 9± is 0.13. The temperature 9± increases with R slightly, and becomes 0.19 
at the boundary, which is sraaller t han the temperature assumed in the Wien fireball model 



( llwamoto fc Takahara II2OO21 . 120041 ). In this case, because of pair-annihilation, the number of 
pairs may decrease seriously before pairs and photons decouple. The temperature 9e is almost 
the same as 9± so that we do not plot 9^. At the same time, the pair temperature does not 
show a adiabatic cooling characteristic of the fireball. This is due to a distributed Coulomb 
heating by protons. If this heating is effective up to the photosphere, pair annihilation 
problem may be significantly alleviated. In this example the proton temperature Qp at the 
center is 260, which corresponds to 1.5 x 10^^ K. On the other hand, at the boundary. Bp 
becomes 1760 (10^^ K). Since this proton temperature is too high to confine by gravitational 
potential of the central black hole, protons will also outflow from the surface, which is in 
contradiction with our basic assumption of a static proton distribution. 

The outflow is due to the anisotropy of the photon field. In Figure 3 we plot densities 
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Fig. 3. — Photon densities of outgoing (solid) and ingoing (dashed) photons for a; = 1 in 
the coordinate frame K. 



of outgoing {n = 1) and ingoing {/i — —1) photons n^{ii,x)dixdx for x = 1 in the frame K. 
In the neighborhood of the center the densities of outgoing and ingoing photons are almost 
the same within errors in our computation. As it goes outside, outgoing photons dominate 
over ingoing photons. This is due to the spherical symmetry in geometry and the beaming 
effect by the mildly relativistic outflow (see also Figure 4). The very low density of ingoing 
photons at the boundary assures the validity of our treatment in computation, in which we 
neglect photons coming from outside of the boundary. Because of the beaming effect, the 
density of outgoing photons is almost constant in the frame although the photon density 
in the comoving frame K' decreases with R. Figure 4 shows the photon spectra of outgoing 
photons for i? = 0, i?o = -Rout/2, and -Rout- The beaming effect is clearly seen in these plots. 



In Figure 5 we plot the radiative force, frictional force, and reaction force due to photon- 
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Fig. 4. — Photon spectra of outgoing photons for R = 0, i?out/2, and i?out in the coordinate 
frame K. 

emission per particle as ratios of timescales; the dynamical timescale ^dyn = Rq/c io the 
acceleration (deceleration) timescale tacc = Tn±mec/ T , where T — F^, Fe + Fp, Frea for the 
radiative, frictional, and reaction forces, respectively. The reaction force causes the reduction 
in the inertia of the flow rather than flow-deceleration. The radiative force also causes both 
the increase in the inertia and flow-acceleration. Thus, Figure 5 does not directly explain the 
growth of U, but wc can see a very powerful role of the radiative force. Around the center, 
the radiative force per particle increases with R because of the anisotropy of the photon 
field. At ~ i?o, tdyn/tacc due to the radiative force attains the maximum ~ 11. Outside 
i?o, the decline of the photon field in the frame K' leads to the dechne of the radiative force 
per particle. In spite of the powerful radiative force, U does not attain the value expected 
simply from ^dyn/^acc ~ 10, because the frictional force competes against the radiative force. 



dyn' ^acc 




Fig. 5. — Radiative (solid), frictional (dashed) and reaction (dotted) forces per particle. 

This strong friction will in turn accelerate the background proton-electron plasma. Although 
we assume that the background plasma is static, in reality it will outflow at a sub-relativistic 
speed. Although U increases monotonically with R, the frictional force decreases outside 
~ Ro owing to the smaller background density and the cross sections of Coulomb scattering 
above [/ ~ 1. 

The energy outflow rate in the form of the pair-fluid is 7.7 x 10^^ erg s~^ in this example. 
The rest of the energy escapes as photons from the boundary. In Figure 6 we plot the energy 
spectrum of the photons escaping from the boundary P^{x)dx. By the analogy from the 
Planck spectrum we can obtain the "photon temperature" from the blue-shifted maximum 
point at X = Xmax = 0.63 in the energy spectrum as 6^ = Xmax/ (2.822r) = 0.15, which is close 
to the pair temperature 6±. Therefore, we could not expect a rise of the pair temperature 
outside the boundary. For reference, we also plot the Wien spectrum oc exp [—x/ {T9^)]dx 
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Fig. 6. — Energy spectrum of the photons escaping from the boundary (sohd). The dotted 
hne is the Wien spectrum for reference. 



in Figure 6. The spectrum has a broader feature than that in the Wien spectrum. Although 
we can see a sign of upscattering of soft photons originating from bremsstrahlungs, there are 
huge amounts of soft photons {x < 10~^) in comparison with the thermal spectrum. 

Despite of the low photon temperature, the positron creation rate (see Figure 7) and 
annihilation rate are almost the same within -Rout- Therefore, we do not plot the positron 
annihilation rate in the figure. Outside the boundary the flow continues to be accelerated 
by the radiative force until the photospheric radius. Although the pair annihilation process 
exceeds the creation process outsid e -Rmit, the number of pairs w ill be frozen at the photo- 
spheric radius. Using the results in llwamoto &: Takahara I (120021 ) . the Lorentz factor at the 

1/3 

photosphere may be T^|^^_ ~ 2, where Tout is the Thomson optical depth at the boundary. 
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Fig. 7. — Positron creation rate. 



The final energy-flux fraction of the outflow to L may be on the order of ~ 10 %. 



For L — 10^^ erg s""^, we have also simulated for another parameter set: Rq — 10^^ cm 
with no = 10^° cm~^. However, the values d±, U, rj, and 9j at the boundary are almost 
the same as those in this example. Therefore, the results may not be so sensitive to the 
compactness of the plasma. 



4. CONCLUSIONS AND DISCUSSION 

This is the first attempt to produce fireballs from hot plasmas, whose size and den- 
sity are characteristic of accretion plasma in AGNs. For microphysics. Coulomb scattering, 
Compton scattering, bremsstrahlung, electron-positron pair annihilation and creation are 
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taken into account. Using the multi-fluid approximation, we obtain spherically symmet- 
ric, steady solutions of radiation and pair outflows for L < 10^^ erg s~^. As the luminosity 
increases, sufficient amount of electron-positron pairs outflow with a mildly relativistic veloc- 
ity. The pair outflows are driven by powerful radiative force rather than their own pressure 
for L > 10^^ erg s~^. Although we had expected that the dynamical effect of the outflow 
would increase the temperature of the pair plasma, in our results of L > 10^^ erg s~^, the 
temperature is rather low to achieve F ^10, di fferently from the assumption in the Wien 



fireball model (jlwamoto fc Takahara 1120021 . |2004| ). However, this difference should not taken 



too firmly at the present stage as is discussed below. 

There are several caveats in our present simulations. First, our calculation regime may 
not be large enough to accurately treat the whole spatial range within the photosphere. 
Distributed heating by protons in the outer region may affect the resultant pair outflow; one 
signature seen in the simulation is the lack of an adiabatic cooling in the flow. Secondly, 
the Coulomb friction rate and the resultant proton temperature suggest that protons will 
also outflow with sub-relativistic velocity. The flow of the background plasma will somewhat 
enhance the outflow of the pair-fluids. If we can raise the temperature to only several times 
the value in our simulation, the fireball can form a relativistic outflow. The temperature, 
density, and velocity of the pair fluid may be strongly influenced by spatial distributions of 
the heating rate and proton density. Qualitatively speaking, shallower proton density distri- 
bution will enhance the outflow, while more centrally concentrated heating rate will enhance 
the photon luminosity. In order to produce a fireball with a sufficiently high temperature, we 
need to continue investigating plasmas with different profiles including the flow of the back- 
ground plasma. While our simulation is one dimensional, two or three dimensional effects 
may be important. 

Another method we should consider is kinematical treatment of the flow. The multi- 
fluid approximation has been adopted in our simulation for simplicity. To carry out more 
accurate simulation, we should solve coUisional Boltzmann equation taking into account pho- 
ton emission. The kinematical treatment may result in a broader distribution of pairs than 
Maxwell-Boltzmann distribution, which raises the effective temperature. The broad photon 
distribution in our simulation suggests such a possibility. The effect of the static electric 
field discussed in §2.3, which may heat the pair-fluid, can be included in the kinematical 
treatment. However, such simulations require enormous cost in computation. In addition, 
photon emission from plasmas with multi-component, which is not in Maxwell-Boltzmann 
distribution, is too complex to evaluate in advance of simulation. Although the kinematical 
methods can include several plasma effects we have neglected, such as two-stream instability 
etc., simulations based on the multi-fluid approximation are irreplaceable way so far. 
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In this paper we have obtained steady solutions for outflows. Since the time scale 
of radiative transfer is much longer than the dynamical time scale, the pair and electron 
temperatures become low. The lifetime of accretion plasma may not be long enough to 
achieve a steady state radiation field obtained here. As the original idea of the Wien fireball 
model, an instantaneous heating may produce a runaway production of relativistic pairs 
because of slow annihilation in a hot plasma. Therefore, time-dependent simulations should 
be examined in a future work. 

The authors appreciate the anonymous referee and M. Kusunose for their helpful advice. 
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